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CN ; Abstract 

"q ' We introduce a generalization of 1 + 1 dimensional large Nc QCD, which we refer to as "reduced" 

■ QCD, or rQCD. In this model gluons and quark momenta live in 1+1 dimensions only, but the quark 

■ spin and all other particles (leptons, and the photon) live in the full 1 + 3 dimensions. The bound 
, states of quarks and antiquarks are identical to those originally described by 't Hooft (except that 

I \ there are new transversely polarized states previously excluded), so the model is exactly soluable. 

\ However, significant differences arise when the model is applied to electromagnetic interactions. 
After reviewing the strongly interacting sector of the theory, we discuss deep inelastic scattering 

■ (DIS) in this model, and show that the new states with transverse polarization give the Callan- 
, Gross relation and remove the pathological features of the original 1+1 dimensional description. 

^ \ We conclude that rQCD gives a satisfactory description of the phenomenology and provides a deep 

O ' understanding of both duality and DIS. 
^ ■ 

Y ! PACS numbers: 12.39.-x, ll.lOSt, 13.40. Gp 
^ . 

c ■ 

> : 

H , I. INTRODUCTION 

In 1974, 't Hooft described the behavior of the large N^. limit of QCD in 1 + 1 dimensions 
P]. He showed that the model confined quarks and was exactly soluable. His work was 
followed quickly by several papers 0, 0, 0] discussing various aspects of the model. The work 
by Einhorn {3] discussed the 't Hooft model predictions for deep inelastic scattering (DIS). 
Einhorn found that the DIS cross section was proportional to the square of the undressed 
quark masses, and hence would approach zero as the quark masses approach zero. This 
pathlogical result made it difficult to apply the 't Hooft model to the phenomenology of DIS; 
clearly someting major was missing. Of course, limiting the physics to 1 + 1 dimensions also 
excludes transverse spin degrees of freedom, so Wi was zero and the Callan-Gross relation 

did not hold. 

Since then QCD in 1 + 1 diminsions has been studied from a variety of points of view 
[li IE and applied recently to the study of duality in heavy meson decays 0, [l2 Jl3l | . 

However, to the best of our knowledge, no one has followed up on Einhorn's study 3| of 
DIS. 
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A few years ago we found that a very simple toy model of scalar and spinor fields in 1 + 1 
dimension gave nice plienomenological results for DIS .14] . By extending the spinor degrees 
of freedom to 1+3 dimensions we were able to recover the Callan-Gross relation. We decided 
to extend our approach to the 't Hooft model, and learned of the pathology in its description 
of DIS. It turns out that this pathology can be removed by extending the spinor degrees of 
freedom to 1+3 dimensions, just as we did before, and this paper grew out of that study. 

This paper is divided into five sections, with seven appendices that include many of the 
details. Following this brief introduction. Sec. II reviews the 't Hooft model results for 
the strong sector. Our discussion focuses on (i) a method of treatment of the confining 
interaction that removes all singularities from the theory and gives finite dressed quark 
masses, and (ii) the consequence of including transverse degrees of freedom for the spin of the 
quarks. We calculate the dressed quark mass, and the properties and spectrum of qq bound 
states, and also study the consequences of the completness of the bound states. We show how 
to construct an off-shell qq scattering matrix, and show that the on-shell scattering matrix 
must be zero, as required for confined particles. Sections III and IV study electromagnetic 
interactions and deep inelastic scatterng. Here the presence of transverse spin degrees of 
freedom completly alters the Einhorn discussion, giving us phenomonologically useful results. 
We conclude with a brief discussion. 

In a subsequent paper [l^ we plan to present a numerical study of the approach to scaling 
in the DIS limit. 



II. THE 'T HOOFT MODEL 



A. The Lagrangian and confining interaction 

We propose that the action have the following form 

<S = Sqcd + Sqed (2.1) 

where the QCD part of the action (gluons and quark momenta or coordinates) live in 1 + 1 
dimensions only, but the quark spin and the QED part (leptons and the photon) live in the 
full 1 + 3 dimensions 

Sqcd = J dtdz CQcoit.z) 

Sqed = j dt dv CQEDit,^) ■ (2.2) 

In his original paper, 't Hooft ^ discussed the strong interactions only, and in this section 
we review the 1+1 dimensional model of QCD that he presented. The QED part of the 
action (j2.2|) . and reduced QCD (rQCD), will be discussed in Sec. IIIII 
The QCD Lagrangian density is 

CqcdH, z) = --Tr [F^^F^.] +Y.q, (iD^Y - ^oO q^ , (2.3) 

i 

where = qi{t,z) is the quark field with flavor i and bare mass moj. The gluon field 
quantities are 

A'^ = Ia^^A. (2.4) 
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F^y = d^,A^ - d^A^ + igo [A^, A^] 
= d^ + igoA^, 

where F^^ is the gluon field tensor and A^ = A'^it, z) are the gluon fields with the Lorentz 
index ^ and the color index a. The gluon fields have components in the and 3 direction 
only, so the sum over the index is restricted to and 3. The QCD coupling constant is go 
and A" are the generators of the SU{Nc) color group, normalized to 



trace [A'^A^ 
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E^°^"= ^ > (2-5) 



where matrix multiplication of the A matrices is implied. As Nc —>■ oo, finite results are 
obtained if go as 

90 = ■ (2.6) 

where the effective coupling (7 is a constant. 

While the space-time coordinates, momenta, and Lorentz vector sums in Eq. ()2.3|) are 
restricted to 1+1 dimension, we assume that the gamma matrices have the usual 4x4 Dirac 
structure. The quark fields are therefore a direct product of four Dirac dimensions and iVg 
color dimensions. 

Note that the action (j2.2p is Lorentz invariant under the subgroup 6*3 of Lorentz trans- 
formations that leave the xy plane invariant. Specifically, the group is generated by the 
hamiltoninian H, and the boost /C3, momentum Ps, and angular momentum jTs, operators 
that generate boosts and translations along the z axis, and rotations about the z axis. Four 
of the six commutation relations between these generators are zero, and the other two close 

[/C3,P3] = -^H [/Ca,?^] = -^^3, (2.7) 

guaranteeing that G3 is a group. 

Following 't Hooft [l| we introduce light cone variables: 

6+ = (6° + }?) (2.8) 

= {bi,h} (2.9) 

for any arbitrary vector b (in the QCD sector, the only perpendicular components come 
from matrix elements of 7_l). Note that the scalar product of any two vectors a and b is 

a^j^b^ = a ■ b = a+b^ + a^b+ — a_L ■ b_L (2.10) 

The derivatives are defined 

Pi — ^ — ^ ( ffi r)A — ^ ( ^ ^ \ 
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so that the divergence of a two- vector is 



dob^ + dsb^ = d-b+ + d+b- . (2.12) 

In the same way, we can define the +, — , and _L components of the 7 matrices. The 
anticommutation relations are all zero except for 

{7+,7-} = 2 

{7x,7x} = {7y,7i/} = -2. (2.13) 

Since the gluon fields are confined to 1+1 dimensions, there is only one nonvanishing 
component of the gluon field strength tensor 

F+_ = = d+A^ - + tgo[A., A+] . (2.14) 

The QCD part of the theory is simplified if we choose the light cone gauge, where A- = 0, 
so that the commutator contained in the field tensor F+_ disappears and there is only one 
nonzero component (A+) of the gluon field. The Lagrangian density ()2.3|) then reduces to 



^QCD = -Tr 



The equation of motion for this field is then 

2 



{d^A+f\ +Yqi (i9+7_ + i9„7+ - ^o7-^+ - "^oi) Qi ■ (2.15) 



d'.AX = 1^-^ AX = -go E 1^ A V Qi ■ (2-16) 

The solution of ()2.16p is 

A+(a;+,x_) =9oJ dy+g{x+-y+) ^ gi(?/+, x_) A''7_ gi(?/+, x_) , (2.17) 

i 

where the Green's function Q is 

'3{x+ - y+) = \x+ - y+\ + ci{x+ - y+) + C2 . (2.18) 

The coefficients ci and C2 cannot be determined without knowing the boundary conditions, 
so they are free parameters. The gauge condition did not ehminate all superfiuous degrees 
of freedom, just as the Coulomb gauge, or the Lorentz gauge, do not determine uniquely 
the photon propagator in QED (Gribov ambiguity). We can therefore set the coefficients ci 
and C2 equal to zero (a specific choice of gauge) in order to simplify our calculations. 

Einhorn discussed the gauge issues related to these two parameters and showed that the 
eigenvalues of the two body bound state equation are independent of the choice of ci and 
C2. However, the dressed quark mass does depend on the choice of C2, and this in turn 
implies that the location of the quark mass pole is gauge dependent. Since any physically 
meaningful quantity is gauge invariant, we can conclude that the location of the mass pole 
is not physically meaningful, and this can happen only if the quarks are confined, so that 
free quark states do not exist. This is an indirect consequence of confinement. Direct 
consequences of confinement will be discussed in the next subsection when the two-body 
bound state equation is discussed. 
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Equation ()2.17p shows that the gluonic field is no longer a dynamical variable, and that 
there are therefore no ghosts. 

The Fourier transform of the Green's function ()2.18p gives us the gluon "propagator" , or 
more precisely the momentum dependence of the effective quark-quark interaction. With ci 
and C2 equal to zero, it can be written 

/oo f 1 fOO rip \ 

^cix+e^(^-^+)^(x+)^|^-5(A:.) /^^j • (2.19) 

This is a singular operator that is well defined only when it is part of an integral over k_ . The 
second term in Eq. (j2.19j) was introduced by Gross and Milana ^] in a different context. 
Its purpose is to preserve the condition 

/oo _ 
dk_g{k^) (2.20) 
-oo 

and it also insures that the potential is finite for any finite value of x+. 

The only variables remaining are the quark fields. The Feynman rule for the (undressed) 
quark propagator is 

- iSo{k) = = -I 2 — KTl — (2-21) 

and the quark- "gluon" coupling is 

- = -^^?o A V , (2.22) 

leading to the following result for the exchange of gluons with momentum k- between quarks 
with Dirac indicies 1 and 2 



= -^9l ^ S i 0127- Qik-) - -t9V-Ou7- Q{k.) , (2.23) 

where we used Eq. fl2.5p (with matrix multiplication of the A matrices implied), and took 
the limit as iV^ ^ cxd using the definition ()2.6p . 

The quark self-energy and the dressed quark mass are calculated in the next subsection. 



B. Dyson-Schwinger Equation for the dressed quark 



We determine the dressed single quark propagator, S{p)^ using the (one body) Dyson 
Schwinger equation (DSE): 



s{p) = So{p)-s{p)j:{p)So{p) 



So{p) - S{p) 



-ig 



(Pk 
(2^ 



Q{p^ - A;_)7_S'(A;)7_ 



So{p) 



(2.24) 



shown graphically in Fig. ^ The gluon interaction does not mix quark fiavors, and the 
calculation is identical for each fiavor of quark, so the flavor index is suppressed. The 
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FIG. 1: Dyson-Schwinger equation for the quark. The corkscrew hne is the gluon interaction, the thin Une 
the undressed quark propagator, and the heavy sohd hne the dressed quark propagator. 



rainbow approximation (undressed vertices and the absence of the quark loops from the 
gluon propagator) is justified in the large Nc limit 0. Since for every internal loop there 
is a factor of = g'^/{2Nc), and a multiplicative factor of I]a A"A" = 2Nc, the color 
dependence disappears. The vertex corrections and the quark-gluon vertices do not have a 
multiplicative factor of Nc, and are therefore supressed in the large Nc limit. 

In Eq. ()2.24p . dPk = dk_dk+, and since D does not depend on it follows immediately 
that the self-energy integral does not depend on either, and must have the form S(p) = 
B{p Hence the dressed propagator is of the form 



Sip) 

where, using 7_7+7_ = 27_, 

B{p^) = -2tg^ 



1 



"^0 - P-1+ - [P+ - B{p^)] 7_ - ie 
the self-energy contribution is then 

dk_dk^ k_ Q{p- — k_) 



(2vr) 



TTln 



2k_{k+ - B{k_)) 



Performing the /c+ integral gives 

dk 



ml-2k_{k+- B{k_)) 
and substituting this back into Eq. ()2.26p gives 



It 



sig{k_ 



dk^Q{p^ — kJ)sig{k. 



Using (E:T!H) for g{p^ - k^ 
B{p-) 



gives 



dk^ 



sig{k-) 



sigip-) 



(p_ — k^y (p_ — k_] 



2np- 



Substituting this back into Eq. ()2.23|) gives 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



Sip) 



mo - p-'j+ 



p+ + 



27rp_ 



mo + P-7+ + \P+ + 



2iTp_ 



7- 



2p+p_ 

TT 



7- 



. (2.30) 
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Note that the quark mass pole has been shifted to a smaller value 



mg — = ml — — , (2-31) 

vr 

where m is the dressed quark mass. It turns out that choosing ci 7^ will not effect the 
dressing of the mass, but choosing C2 7^ would give 



= ml - — {1 - TTC2\p-\) . (2.32) 



TT 



This gives a mass that is momentum dependent and not covariant. Hence the dressed mass 
is gauge dependent and unphysical. With these cautionary remarks, we choose ci = C2 = 
because the mass is covariant. One undesirable feature of this choice is that < in the 
chiral limit. This could be corrected by choosing C2 large and positive. 

Having obtained the dressed propagator, we are able to proceed with the two body bound 
state calculation. 



C. Two-body bound states 

1. Algebraic form of the two body equations 

Consider a bound state of a qq pair. The quark has dressed mass mi and electric charge ei, 
and the antiquark (which might be of a different flavor) has dressed mass m2 and charge 62- 
The momentum of the bound state is r, the momentum of the quark is p and the momentum 
of the outgoing antiquark is r—p, as shown in Fig.|21 In much of the following discussion, we 
will treat the outgoing antiquark as an incoming quark with momentum p — r, and following 
this convention label the bound state vertex function r{p,p — r). For a color neutral state 
we may carry out the color sums using (j2.23p . giving the Bethe-Salpeter equation for the 
bound state vertex function: 

T{p,p-r) = I J 0-^V{p,k)^^S^{k)T{k,k-r)S2{k-r)^_, (2.33) 

where 5*1 and 5*2 are dressed quark propagators, and the kernel is the singular operator 

V{p,k)=g^g{p^-k^). (2.34) 

With the substitution T{p,p — r) = 7_ G{p,p — r) [1], Eq. ()2.33p becomes: 

^ _ , . /• <Pk V{p, k) k_ (fc_ - r_) G{k, k-r) 
U[p,p r)-^i] ^2^^^ _ ^2 _ _ _ _ . (^.3bj 

The formula shows that G{p,p — r) does not depend on and hence the k_^. integration 
can be carried out immediately. 

Integrations of the variable k^ over two quark propagators occurr frequently, and are 
worked out in detail in Appendix ^ It is convenient to introduce the momentum fractions 

y = ^ z = ^. (2.36) 
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FIG. 2: The Bethe-Salpeter equation 



If the fraction z lies in the interval [0,1], Eq. (j2.35|) reduces to 

'tti + l a2 + l\ ^ , , /-i $„(?/, r) - $„,(2;,r) 



[ii^n{z,r) = + - <^n{z,r)-f dy- 

\ z 1 — z J Jo 



{z - yY 

= -^E{z)^^{z,r) (if ^€[0,1]), (2.37) 
9 

where / is the principal value integral, and ai, and 02 are dimensionless parameters 



ttM^ vrm? -ami , ^ 

/^n = ^ > «i = ^ ' «2 = ^ , (2.38) 
9 9 9 



with = the anticipated mass eigenvalues that will emerge from the solution of the 
equation, and rrii the dressed quark masses. The two-body "wave function" $„(2;, r), is 
defined by 

ATr, Gn{z,r) ^ Mr_ Gn{z,r) _ Mr_ Gn{z,r) 

\ Z I - Z J TC \ Z 

where Gn{z,r) = G{p,p — r) is the anticipated eigenfunction solution for the vertex cor- 
responding to the nth bound state, and jV is a normalization constant. It is chosen so 
that 

'dz^l{z,r), (2.40) 





and will be calculated in subsection IIII CI Note that Eq. ()2.37|) defines the two-quark 
Hamiltonian, H{z), on the interval < z < 1. The bare quark masses enter into the 
Hamiltonian, while the dressed quark masses enter into the relation ()2.39|) between the 
vertex function and the wave function. 

If the momentum fraction z lies outside of the interval < ^ < 1, the vertex function 
defined by Eq. ()2.35|1 is not zero, but can be obtained by quadrature from the vertex function 
defined inside of the interval [0,1]. Using ()A8|) and the definition (|2.39p gives 

z ^ 1-z ^" 
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These definitions and equations will be used in the following sections when the hadronic 
currents are computed. 

An alternative form of Eq. ()2.37p follows if we integrate the last term by parts: 




z-y) 

l,r) - ^n{z,r)' 
I- z 



(2.42) 



where $^(?/,r) = d^{y,r)/dy. If $„(0,r) = $„(l,r) = 0, which is true in all but the chiral 
limit, this form of the equation follows directly from 't Hooft's principle value prescription 
fT], so that our two-body equation is identical with 't Hooft's, even though we have finite 
dressed quark masses. 

To avoid confusion later, we call attention to the fact that the wave function defined by 
Eqs. (|2.39p and (|2.41|) is not identical to the Bethe-Salpeter wave function. The latter can 
be defined 

n,,r)^ ^ ,^^-'P-;-'-/'f-P-;> .. . (2.43) 

[mf — — te) [1712 ~ [p ~ '^r ~ "^^J 

where the factor p-{p — r)_ in the numerator comes from the 7_ factors in the quark 
propagators S. This wave function carries the p+ dependence of the quark propagators. If 
we integrate (|2.43|) over p+ (using the fact that G does not depend on p+) we obtain 

Hp-,r)^ dp+^ip,r) = \ A{z,r^) "^^^'^^ it ^ G [U, ij _ ^2.44) 
■^-^ [O \izi [0, 1] 

This is the "equal wave function, and equals the wave function we are using only in 
the region z G [0,1]. As we have seen, the wave function <l>(z, r) that solves the 't Hooft 
equation is non zero jor all z. 

We conclude this discussion of the two body wave equation by returning to Eq. ()2.33|) . In 
place of the replacement r(p, p — r) = 7_ G'„(2;,r), it could equally well have been assumed 
that 

r(j9,j9 - r) = 7_7fcG„(z,r) , (2.45) 

where k is either x or y (recall that we are assuming that the Dirac matrices live in 1+3 
dimensional space). Note that the Bethe-Salpeter equations with the substitutions (|2.45|) re- 
duce to the same equation (j2.35p for and hence each two-body state is triplely degererate, 
with spin structure given by 7_, 7_7x, or 7-7y. In each of these three cases the wave func- 
tions and bound state mass are identical. The additional states with the jx or 7^ structure 
are outside the scope of the original 't Hooft model, but must be considered in reduced 
QCD. They are states with quark spins in the transverse direction, and will play a major 
role in the discussion of DIS below. 

The chiral limit is exactly soluable and of some interest. If = 0, then Eq. ()2.37|) has 
the normalized ground state solution 

Mo = 0, %{z,r) = l. (2.46) 

Outside of the interval [0, 1], Eq. ()2.41|) gives the same result, so that 

Mz,r) = l (2.47) 

for all z in {—00, 00). 
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FIG. 3: Diagrammatic representation of confinement via two-body interaction. The x indicates that the 
dressed quark is on-shell, and at the point when both quarks are on-shell, the vertex function is zero. 



2. Confinement in the 't Hooft model 

The 't Hooft model illustrates how confinement can be realized in two (apparently dif- 
ferent) ways. If the singularities of the confining interaction are left unregularized, as they 
were in 't Hooft's original paper [H, the quark masses are infinite. Even if the singularity is 
regularized as in this paper, whenever < g'^/ir the dressed quark mass will be unphysical, 
and there is no real quark mass pole. This insures that the quarks are confined. 

In the second case {ml^ > j-a and > 0) the dressed quark mass is real and positive 
and the there is a quark mass pole. What keeps the quarks confined in this case? The 
answer is that the quark vertex function Gn(^, t) has a zero at precisely the value of z where 
all of the quarks in the bound state could be on mass shell. This requirement is illustrated 
diagramatically in Fig. El and seems miraculous. However, as discussed in Ref . , it is 
only a consequence of fact that the wave junction is finite at the pole (because the wave 
equation permits no singularities), and therefore the vertex function 



must be zero. Precisely the same mechanism works in the Schrodinger description of two 
nonrelativistic particles bound by a linear (or any other) confining potential. 

In summary, quark mass poles are always unphysical because the dressed mass is gauge 
dependent, and bound state poles are physical because they are gauge invariant. In this 
respect the undressed quark mass is the only physically significant parameter. The two 
ways of viewing confinement are therefore not exclusive, but equivalent, since they are 
related through gauge transformations. 

3. Numerical solutions for meson wave functions 

The bound state equation was solved numerically using modified cubic splines. The 
splines and the technique are fully described in Appendix [HI Each spline is a smooth 
function (continuous up to its third derivative) and has support over a small region of the z 
axis. The splines overlap, so that by independently adjusting the height of each spline it is 
possible to smoothly approximate any function defined on the interval [0, 1]. 

We have modified the standard spline technique to insure that the numerical solutions 
satisfy the boundary conditions exactly. These boundary conditions are determined by 



examing Eq. (j2.37j) in the vicinity of z ~ and z ~ 1 lll|, and require the wave function 
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TABLE I: Numerical values of for moi = rno2 = 1-5 and various (number of splines) and Ug 
(number of guass points for each spline segment). 

level Us 20 40 60 

n Ug 2 A 6 246 246 

1 13.730 13.732 13.732 13.717 13.718 13.718 13.714 13.715 13.715 

2 25.613 25.614 25.615 25.572 25.573 25.573 25.563 25.564 25.564 
10 108.95 109.35 109.36 108.42 108.43 108.43 108.28 108.28 108.28 
20 186.17 210.60 210.58 208.82 209.69 209.69 208.88 208.97 208.97 
30 _ _ _ 299.64 314.02 314.02 308.04 309.38 309.38 
40 _ _ _ 362.04 409.77 409.74 410.31 411.24 411.24 



TABLE IL Numerical values of for m-oi = mQ2 = 0.5 and various and Ug. 
level Us 20 40 60 
reg 4 6 8 4 6 8 4 6 8 

1 2.6196 2.6625 2.6841 2.6485 2.6786 2.6936 2.6616 2.6858 2.6979 

2 10.540 10.650 10.706 10.616 10.692 10.731 10.648 10.710 10.741 
10 86.170 86.371 86.477 86.055 86.288 86.409 86.197 86.400 86.503 
20 199.44 199.43 199.43 183.93 184.11 184.21 183.53 183.77 183.90 
30 _ _ _ 287.16 287.22 287.25 282.09 282.27 282.36 
40 _ _ _ 399.18 399.15 399.15 333.13 333.22 383.27 



go like 

$„(z,r) z^^ (as z ^ 0) with 7r/3i = —ai tan(7r/3i) 

<l>„(z, r) ^ (1 - zf'' (as z^l) with vr/^a = -«2 tan (77/^2) . (2.49) 

Note that the = in the chiral hmit (when ai = —1), consistent with the exact solution 
()2.47p . To fit these boundary conditions the splines closest to z = and z = 1 have these 
fractional powers built into their functional form, as described in Appendix O 

The wave function is expanded into Ug splines, each with four segments (except for those 
at the boundary, which have only three). The splines overlap, and taking into account that 
there must be at least three splines (one regular spline in the center and one modifed spline 
at each end), the total number of segments is + 1, dividing the interval into segments of 
length l/{ns + 1). The equation is turned into a matrix equation by integrating over each 
segment using ng gaussian points. Tables HI and ITTl show the numerical stability. If the quark 
mass is larger than unity, better than 1% accuracy is achieved with only two gauss points 
per interval, provided n ^ Ug. As the level number n increases above ^ ns/2 toward n^, the 
reliability of the calculation decreases. To study highly excited states, it is sufficient to use 
Ug = 2, and push as high as possible. 

The situation is somewhat different if the bare quark masses are less than unity (Table 
EH). In this case the wave function is very steep at the boundaries and 4 or 6 gaussian 
points per interval are needed to get 1% accuracy for the low lying states. In a subsequent 
paper (l^, we will limit our numerical discussion of duality to cases with moi > 1 where the 
solutions are very stable, and there are no poles in the region outside of [0,1] [see the brief 
remarks following Eq. ()F14|1 . 
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FIG. 4: The square of the bound state mass (in units of g^/7r) versus the order of the state n, computed 
using 60 spHnes. The straight hne is the hnear relation mrg^ . Note that the masses follow this linear relation 
up to near n ^ Ug = 60, where the departure from linearity is due to numerical inaccuracy. 
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FIG. 5: The ground state wave functions for equal bare quark masses mo of 1.1, 1.5, and 3 (the wave 
functions become steeper at z = and 2: = 1 as the quark mass decreases). In all cases 60 splines were used, 
and the wave function was normalized to a maximum value of unity. 

Figure m shows the linearly rising Regge trajectory characteristic of the 't Hooft modeL 
Figures IHE show the ground, first excited and 10th level for three cases with equal quark 
masses of different values. Note the "tails" of the wave functions outside of the region [0,1], 
given by Eq. (j2.4H) . These tails are large for small quark masses, but shrink quickly as 
the quark mass increases, and are also small for highly excited states. Figure |H1 shows the 
20th level for spline numbers 20, 40, and 60. Note that the state is not fully described even 
with 40 splines, supporting the observation that an accurate description of states of order n 
requires a spline number of approximately ~ 2n. 
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FIG. 6: The second state (first excited state) for the same three quark masses shown in Fig.O (See caption 
to Fig. El) 








z 

FIG. 7: The tenth state for the same three quark masses shown in Fig.|31 (See caption to Fig.|31) In this 
figure the dots are sin(107rz), showing that this is an excellent approximation for small quark masses. 

D. Scattering of quarks and antiquarks in the presence of confinement 

1. Completeness relation and the two-body Greens function 

The Hamiltonian defined in Eq. ()2.37|) is hermitian on the interval [0, 1] and it is straight- 
forward to show that the different solutions are orthogonal on this interval. With the nor- 
malization given in Eq. (j2.4Uj) 

/ dz^n{z,rn)^m{z,r^) = 6nm. (2.50) 

Jo 
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z 1 

FIG. 8: The 20th state for and equal quark mass of 1.5. The small amplitude curve was calculated for 
Us — 20, and completely fails to represent the state. The dotted and full amplitude solid lines are calculated 
for Us = 40 and 60, respectively. Note that — 40 does well except for the first and last oscillation. 



The completeness of the eigenfunctions imphes that 

"^niz, r„)$„(2', r„) = 6{z - z') if z, z G [0, 1] . 



(2.51) 



The Greens function is therefore 



G{z,z\q^) = Y. 



(2.52) 



where q is any external momentum. The completness relation shows that G has the property 
of a Greens function 



{q^-H{z)\G{z,z\q^) = E " ^ 



$„(z,r„)$„(z',r„ 



g2 - Ml 
= 6{z - z') if z, z' G [0, 1] . 

The same argument also gives 

G{z, z\ g2) {g2 _ E{z')\ = 6{z - z') if z, z' G [0, 1] 



[2.53) 



(2.54) 



2. The scattering matrix 

In preparation for the study of electromagnetic interactions, it is useful to introduce the 
qq scattering matrix, defined by the following infinite series 

M(p', p; q) 7i_ ® 72- = V{p\ p) 71 (g) 72- 

+ ^ J 0^,V{p\k)V{k,p) [7_5i(A:)7-]i®[7-52(A:-g)7-]2 
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+ ••-. (2.55) 

The momentum are defined as for tlie two-body bound state, except tliat liere q is tlie 
unconstrained momentum of tlie qq pair and q^ is not an eigenvalue. The kernel V is the 
singular operator (|2.34p . The singularities of this operator will never be exposed because we 
will limit application of this scattering equation to cases where we integrate over either the 
initial or final momentum. Note that the structure of the equation anticipates the reduction 
of the dressed propagators 



so that the factor of 7i_ ® 72- is common, and the series ()2.55|1 reduces to 

d^k V{p',k)V{k,p) 



M{p',p;q) = V{p',p) + ^ j 



d2{k_, k+,q) 



^ ^2 ifK M^fc V{p',k')V{k',k)V{k,p) ^ ^2 57) 

TT^ i TT^ d2{k'_,k'+,q) d2{k-.,k+,q) 



where d2 was defined in Eq. flA3|) . From this result we can conclude that M does not depend 
on p'_^ and p+. 

In applications below, we will encounter this scattering series in the following form 



{MGO) ip'; q) 



d^pM{p',p]q)0{p-,q) f d^pV{p' ,p) 0{p-,q) 



+1 



7r2 d2(P-,P+,g) J 71-2 d2{p-,p+,q) 

d'^k V{p',k) f d^pV{k,p)0{p.,q) 



n"^ d2{k_,k+,q) J tt^ d2{p-,p+,q) 



^ rd'k' rd'k V{p\k')V{k',k) rd'p Vik,p)Oip.,q) 
J Tc"^ J 71^ d2{k'_,k'j^,q) d2{k_,k+,q) J ir^ d2{p-,p+,q) 

where 0{p-, q) is an operator that, by assumption, does not depend on p+. In Appendix O 
we show how to use the Greens function to write this series as 

/' dz M{z\ z- q') = E A(/, M„2) %^ /' dz <|.„(z, r„) 0(^, q') , (2.59) 

Jo L\[z,q ) „ iw^ q Jo 

where we replaced {q-Y M {p' , p; q) / n — > M{z',z;q'^) and 0{p-,q) 0{z,q'^). 

This equation has a nice physical interpretation. It displays the "scattering amplitude" 
as the sum over the propagation of the bound states, which is all confinement will allow 
(without meson decay mechanisms, which are ignored in this paper, no cuts are possible). 
The only singularities in that can occur are poles at the bound state masses M^. 

With these tools in place we are able to study the electromagnetic interaction of hadrons. 
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III. REDUCED QCD AND THE ELECTROMAGNETIC SECTOR 



We return to the discussion of the Lagrangian density introduced in Eq. fl2.2|) and look at 
the electromagnetic sector. In rQCD the electromagnetic fields (and the electron currents) 
are extended to 1+3 dimensions. The Lagrangian density is therefore 

CQEDit,r) = -^F^''{t,r)Fof,^{t,r) + "^e^ qi{t, z)'y^qi{t, z) A^{t,r) 

i 

- ej;(t,rX(t,r), (3.1) 

where y4o(t,r) is the electromagnetic field (with the subscript distinguishing it from the 
gluon field), the vector sums are now over 4 dimensions with yU = {0,1,2,3}, Cj is the 
electromagnetic charge of quarks with fiavor z, J^(t, r) is the electron current that produces 
the electromagnetic field, and 

F^"" = d^'A^o - d^A!^ , (3.2) 

is the electromagnetic field tensor. Because the electron current exists in four dimensions, 
all four components of the electromagnetic field will be non-zero, in general, and all four 
components of the quark current will be excited by electromagnetic scattering. 

We now apply this discussion to electron scattering, where = —q^ > 0, so that 
g_ < 0. The initial particle (quark or hadron) has momentum p = p' — q and the final 
particle momentum p', with momentum fractions 

z' = pL = / _ 1 . (3.3) 

q- Q- 

If the final particle is on shell, then p'_ > and z' < 0, and z' — > — oo as 0. We will 

use the notation jg to denote the bare quark current operator, which is 

Jo^ = Y . (3.4) 

While the dressing of the strong quark-gluon vertex is supressed by the large Nc limit 
P], the quark electromagnetic vertex can be dressed by gluon exchanges. The longitudinal 
part of the dressed electromagnetic vertex {j^) can be computed directly from the Ward- 
Takahashi identity 

q,j''{p',p) = S-\p)-S-'\p'), (3.5) 
where p + q = p'. Using the dressed propagator ()2.30p . we find that 

j- = 7- 

3+= 7+ - 7- 7^ = 7+ + 7- 7^ 7\ ■ (^.6) 

2n p_p- 2n q_ z'[l — z') 

Consequently the j_ component is unmodified. 



A. Electromagnetic coupling to quarks and the quark form factor 

It is instructive to obtain the full quark current directly from the dressing of the quark- 
photon vertex. The vertex is dressed by successive gluon exchanges, as shown diagramati- 
cally in Fig. El When coupling to a photon, there is no fiavor change, and hence, for each 
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FIG. 9: Diagrammatic picture of the gluonic dressing of the quark current. 



flavor i, this series is 

jiip',p' -q) = Y + i 



(Pk 

(2^ 



+ ■ • 



(fk 



(27r)2y (27r) 



2 V{P\P) V{.P, k) 7_ S^ij)) 7_ S^{k) Y Si{k - q) 7_ Siij) - q) 7_ 



(3.7) 



This series is evaluated in Appendix using the methods developed in Appendix O The 
final results are conveniently expressed in term of a quark form factor, defined to be 



M2 + g2 J, 



(3.8) 



where z' = p'_/q- and = —q^ in anticipatation of applications to electron scattering. 
Note that the form factor is expressed as a sum of bound state poles at = —M^, arising 
from the propagation of the bound states that couple to the photon. Vector dominance is a 
rigorous consequence of this model. 

Using this form factor, the currents reduce to 



ji+ip',q) = 7+-7- 



jT(^',9') 



1x - Ix 7- 



g 1 ^ <i 

f Q- 9 



,2 1 



+ ^F,{z',Q'' 



niQi 27[ p'_p- mo 



(3.9) 
(3.10) 

(3.11) 



Note that the terms independent of the quark form factor F^ are identical to the expected 
result ()3.fij) . and that the contribution from the form factor is purely transverse, in that 



q+ji-{p,q) + q-ji+{p,Q) 



0. 



(3.12) 



This result mimics the method for insuring current conservation developed in Ref. [18[, and 
gives some evidence that that method is dynamically sound. 

These results are used in the next subsection to calculate the hadronic transition current. 
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FIG. 10: Currents for the electromagnetic transition from the initial state i to the final state /. The left 
hand diagram is the contribution from the electromagnetic interaction with the quark with mass mi and 
charge ei and the right hand diagram is the same for the antiquark with mass TO2 and charge 62- 



B. The transition current 



In the presence of confinement, and without any meson decay mechanisms, the quarks in 
a qq bound state must remain bound, even after they absorb an energetic photon. In this 
section we calculate the transition current for the process 7* + Mi — ^ Mj, where Mf is the 
mass of a (possibly excited) final bound state and Mj is the mass of the initial (ground) 
state. 

The transition current consists of the sum of two contributions: one in which the photon is 
absorbed by the quark with dressed mass mi, and a second in which the photon is absorbed 
by the antiquark with mass m2- The momentum used to label each of these processes are 
shown in Fig. the initial bound state has momentum P and the final state has momentum 
Pf = P + q. The transition current coming from the first term is 



{f\J'^iPf,P)\^) 



eitcrm 



xTr 



ieiNc 



7- 



1 



k 7x J 



(27r)2 

S^iPf + k)j'^S,iP + kh^S2ik) 
{-Pf - k, -k) Gi (P + k, k) 



(2vr) 



{Pf + kf - ie] [ml -{P + ky - ie] [m: 



(3.13) 



where is the quark current operator for quark 1, as worked out in subsection IIII A[ A^'^ 
is the spin dependent numerator discussed below, and we have allowed for the possibility 
that the final state may have either the 7_ or the 7_7a; structure discussed above (we will 
not consider the 7-7j/ states; they are needed only if we consider the y component of the 
transverse current and give results identical to the results for the x components). Note the 
factor of Nc coming from the sum over all possible quark colors. 

Only the term in the numerator of the S2 propagator equal to /c_7+ can give a nonzero 
result, but the terms that contribute from the Si propagators depend on the matrix form of 
the final state, and of . Because the trace of a product of gamma matrices with only one 
factor of 7x or jy is zero, the transverse currents (with or 7^) will be zero unless the final 
state also includes (or '-fy for the y component of the current). Hence we have the rules: 
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(i) x-type transverse currents couple 7_ states to 7-7a. states only, and (ii) longitudinal 
currents couple 7_ states to each other. Because of the exact degeneracy of these states, the 
results in both of these cases use the same momentum-space wave functions. 

The integral (|3.13|) is evaluated in Appendix |F| The result for the — component of the 
current is 



ei term 



IGeiAT, P_ 



X 



where 



q- 

y = -p: 



+ 1 njd^ ^j{e,Pf)M^,P) 

l + F,{r],Q')] , 



V 



y 



(3.14) 



(3.15) 



[for definitions of all the momentum fractions see Eqs. ()F1|) and ()F2|) ]. the state /_ has the 
same structure as the initial state (assumed to be 7_), Fi is the quark form factor defined in 
Eq. ()3.8p and TZ is given in Eq. ()F14|) . The plus and transverse components of the current 
are 



ei term 



{f.\J+{Pf,P)\i.) 



ei term 



(27r)Ar' 



WeiNcP- 



q+ [mi-v{l-r])Q^F^ir],Q')] 
moi ?7 (1 - r/) 

g+ [ml-r^{l-r^)Q^F,{r^,Q^)] 
q_ T] {1 - rj) Q"^ 



(3.16) 



where the state fx has the structure 7- 7a;. These results are used in the next subsection, 
and in the discussion of DIS. 

Using the definitions of the momenta shown in Fig. ^[ it is easy to see that the second 
term in the transition current can be obtained from the first by letting 1 ^ 2, as discussed 
in Appendix |F| The exact expression for the x component of the transition current, which 
is the sum of the ei and 62 contributions, can therefore be written 



where the reduced current Jjf{Q'^,y) is 



Jir{Q\y) 



i-y 



+ I nAd^ 



with {6, ei, m} = {e, C^} and {6, ^2, 772} = {1 - 1 - e', 1 - r/} and 

1 



Oi 



mi-r]{l-r])Q'F,{r],Q^ 



(3.17) 



(3.18) 



(3.19) 



Finally, the function TZi = 7^(mi, 777,2), was defined in Eq. ()F14|) : the substitutions that 
convert the ei term into the 62 term therefore give 7^2 = 'R.{m2,mi). 
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C. Elastic scattering and normalization of the wave function 



If the scattering is elastic, x = 1 and the momentum fraction y = —q^/P^ must be 
computed from the exact expressions for q and P in the cm. frame 1^. The result, for 
elastic scattering from the state n, is 



y 



+ QJAM^ + 



(3.20) 



As — s> 0, ?/ — » and i] [defined in Eq. fj3.15|) ] approaches — oo, FiijjjQ"^) — > 0, and the 
exact result for the elastic curent becomes 



{n\J^{P,P)\n) 



/ rfe[<^'n(e,^)]' = 2eiP_^„(0), (3.21) 



ei term {2ts:)M 

where charge normalization requires that the form factor J^niQ"^) be unity at = 0. Hence 



AT 




(3.22) 



gives the normalization condition ()2.40|1 . 

Furthermore, elastic scattering requires that, as 0, 



(r/-l)g_ = eP- 



(3.23) 



and hence 



{n\J+iP,P)\n) 



TTl [' 

= 2eiP+--| / di 

eiterm M„ JO 



(3.24) 



= 2eiP+J-„(0) =2eiP+, 

which follows from the wave function identity ()B4|) . Note that the charge is properly normal- 
ized and conserved for both components because we used the correct dressed quark current 

Finally, we point out that the elastic transverse current is zero for all Q^, as required by 
the summetries of the states 



{n\J,iPf,P)\n) = 0. 
We now turn to a discussion of deep inelastic scattering. 



(3.25) 



IV. DEEP INELASTIC SCATTERING 



A. Extracting the structure functions from the DIS cross section 



Using the standard conventions (also defined in Ref. [l^); the DIS cross section is 



_fa_ 
dVt'dE' 



1^2 cos^ (^+2W,sm' 



W2 + 2Wi tan^ 



(4.1) 
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where the structure functions Wi and W2 are part of the hadronic tensor, defined by 




where M = Mi is the target, or ground state mass. The structure function Wi is immediately 
extracted from the xx (or yy) component of the tensor 

MW, = ^^Fiix) = ^f{x). (4.3) 

Here we anticipate that MWi scales in the deep inelastic limit to the function Fi{x) of the 
variable x = Q"^ /{2P ■ g). Note that Wi must be identically zero if the x (and y) components 
of the current are zero. 

The structure function W2 can be extracted from the + and — components of the current 
using 

= \ (W^oo - W^..) = 2nM 1 (1 + -W^Y (4.4) 

In the deep inelastic limit this gives 

UW2 F2{x) = 2xFi{x) + -xW_+ . (4.5) 

TT 

The contribution from W-+ will be shown to go like in the following sections. Hence, if 
Fi is non-zero, we obtain the Callan-Gross relation. However, if the transverse components 
of the current are omitted, Fi = and h'W2 does not scale. In this case it is the quanity 
i'^W2 that scales, and we recover the results of Einhorn j^. 

We now turn to a calculation of the structure functions. We calculate them first (i) 
in the partonic picture, where all final state interations are ignored and the quarks in the 
final state are assumed to be free, and then (ii) in the hadronic picture, where confinement 
ensures that the only possible final states are the bound gg states we have already discussed. 
We emphasize that the partonic picture can never actually occur because the confining 
interaction, which also acts in the final state, can never be ignored. 



B. DIS in the partonic picture 

Consider deep inelastic scattering (DIS) in the partonic picture, which assumes a final 
state composed of free quarks with no interaction. This picture is illustrated in Fig. ^2 The 
full result is the sum of a contributions in which the photon is absorbed by the quark with 
charge ci, and one in which it is absorbed by the antiquark with charge 62- In the center of 
mass of the ejectiles, either particle can go forwards or backwards and the unpolarized cross 
section is the sum over both helicities of the particles. 

Begin by calculating the transverse component (chosen be be x) of the current. The 
outgoing quark has a mass of mi and momentum pi in the center of mass frame of the 
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— 




Pi 

FIG. 11: One of two diagrams that contribute to DIS in the partonic picture. Note the labchng of the 
momenta. 



ejectiles, while the antiquark has a mass of m2 and momentum of p2 (in this section p2 will 
be the physical momentum of the outgoing antiquark, opposite in sign to that used in the 
earlier sections of this paper). The current is the sum of two diagrams, only one of which is 
shown in Fig. ^2 Since the qq vertex contains the factor of 7_, only the factor m +p_7+ in 
the numerator of the dressed propagator ()2.30|) survives, but we will retain only the factor 
of p_7+ which dominates the result at large Q^. The following result is obtained for the 
X component of the current for the disintegration of the ground state (the subscript n = 1 
usually attached to the ground state will be omitted in this section): 



4 = -ei [M(pi)7^7+7_t;(p2)] 
+62 [u(pi) 7-7+7,; t;(p2)] 



jpi- q)-G{pi- q, -P2) 

mf — {pi — qY 
{P2- q)-G{pi,-p2 + q) 



^2 - {P2 - qY 

Introducing the momentum of the bound state, P = pi + P2 — q, and defining 

(Pi-q)- {P-P2)- 



(4.6) 



(4.7) 



[as in Eq. ()2.36|) ]. the first denominator can be written 



^1 - (pi - of 



mi 



{P-P2f 



mi 



M2 

2P_ 



mo 



2p2 



{P-P2)- 



rrin 



1 - z 



2 2 
ml mi 

z 1 — z 



(4.8) 



Here we have made use of the constraints = 2P_P^ and m\ = 2^2-^2+ • A similar result 
follows for the second denominator if we define a new momentum fraction z' 



{P2 -q)- ^ {P -Vi)- 



(4.9) 



With these definitions we obtain 



Jx = -ei[u{pi)-f^-f+-f.v{p2)] 



P-G{z,P) 
A(^,M2) 
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r / N / M P-G{1-Z',P) 

+e2[M(pi)7-7+7xt^(P2)] A(l-z^M2) 

+ ^ [wbi) 7-7+7. ^^(^2)] $(1 - P) , (4.10) 

where the answer has been expressed in terms of the wave function (j2.39|) . 

The current and the DIS cross section are evaluated in the cm. frame, using the kine- 
matics defined in Appendix |E| The unpolarized cross section is obtained by squaring the 
current and summing over spins using 

J2 HPihxl+l-v{p2) v{p2h-l+lxu{pi) = pi+p2-trace [7+7_7+7^7_7^7+7_] 

spins 

= 32pi+p2-, (4.11) 

and neglecting the interference term (which vanishes in the large limit). Observing that 
z z' X in the DIS limit, and summing over all colors of the outgoing quarks (which do 
not interfere) gives 

= ^j^{eipnp^i^'i^,p)+eip?ipn^'{i-x,p)} 

= AnQ'[^) {el ^'{x, P) + ej $^(1 - x, P)} , (4.12) 



spins 

'1 — X 



where the momentum are evaluated in Appendix|El The hadronic tensor defined in Eq. ()4.2p 
is related to the square of the currents by 



W^u _ r dpidp2 
AttM ~ J (27r)2 4pioP2o' 



The Wxx component of the tensor is 3| is therefore 

4b.|(pl0+P20) stii^s 

= 27T {el ^\x, P) + el $2(1 - X, P)} . (4.14) 
This gives the familiar parton model result for the structure function Wi 

MW,= Pi(x) = ^ {el <!>\x, P) + el <!>\l -x,P)} = ^f{x) . (4.15) 

The structure function W2, given in Eq. (j4.5p . is then 

uW2 = 2MxWi = 2xFi{x) = F2{x) = xf{x) , (4.16) 

where we anticipate the result W^^ as 00. This is the Callan-Gross relation. 

Now look at the plus and minus components of the current. If the Dirac space were 
restricted to two dimensions, only these components of the current would exist. The minus 
component can be obtained immediately from ()4.10j) by replacing 7^ with 7_, giving 

J_ = - [uip^) 7- vip2)] ^ {ei$(^, P) - e2$(l - z', P)} , (4.17) 
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The plus component could include the two contributions from the quark current, Eq. ()3.6p . 
but in the partonic picture the outgoing quark is not dressed, so the second part of the 
current proportional to g"^ should be ignored. However, in order to study the contribution 
of this term we will retain it for now. The leading contribution from both terms in the J+ 
current is 



J4 



moi— [u{pi) 7+ 7_ v{p2)\ + 



+ 



62 



mo2 — [m(pi)7-7+^(P2)J + 



TT (1)/ (1) 

2 ^iz,P) 



(2) 
P2- 



P2-{P2- - Q- 



Hpi) 7- v{p2)] 
[u{pi) 7_ v{p2)] 



(4.18) 



The W. 



component of the hadron tensor is then 
1 



W. 



4 \Pz\iPio + P20) 



X 



spins 



(4.19) 



spins 



Using ()4.17p and ()4.18p . and summing over colors, the ef contribution to the current sum 
becomes, at large Q^, 



spins 



32 el Nc 
A/'2 



<^\x,P) 



1 — a; 



P2- 



X 



(1) 
Pi- - 

<l>2(x,P) 



01 



Hence 



lei term 



Ipim? $^(^,P)=4vr^/(x)U 



term ' 



(4.20) 



(4.21) 



and this term is subleading at high . Note that the effect of neglecting the term pro- 
portional to in the + component of the quark current (as we are instructed to do in the 
partonic picture) is to replace m\ by the bare quark mass m^^. If W\ = 0, then i''^W2 scales 
as 



u'^W2 



ei term 



-w. 



^<1 



/(^)Li 



term ' 



(4.22) 



27rM --'^it^™ M 

Up to a factor, this is the result originally obtained by Einhorn Q. In that paper, the 
Dirac space was restricted to 2 dimensions, Wi was necessarily zero, and i''^W2 scaled as 
(I4.22|) . This distribution amplitude does not have the correct physical interpretation. In 
particular, it predicts that the DIS cross section depends on the square of bare quark mass, 
a clearly unphysical result. Only by considering the full four dimensional Dirac space are 
we able to obtain scaling for vW2-, and the familiar physical result ()4.15|) for f{x). 

Another consequence of the fact that is subleading is that we can avoid coming 

to grips with the fact that the current components J„ and J+, as defined above, do not 
conserve current. Since these components do not contribute to the final result, we need not 
discuss how they can be redefined in order to conserve current. 
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C. DIS in the hadronic picture 



1. Exact results for inelastic scattering 

The exact hadronic tensor (j4.2p for electroproduction of a single final state / with mass 
M„ (denoted by n) from the ground state i with mass Mi = M (denoted by 0) is 

/dP 
(27r)2P/o spin, 

= <5(M„ - Po - go) ^ ^ (0|J-Jn) (n|J,|0) . (4.23) 



spins 



This is a deha function. In any physical measurement, the detectors will accept a finite 
range of values of the final electron energy, E' , requiring that this theoretical cross section 
be averaged over E' . This averaging process will be discussed in subsection IIVDI below, 
where the deta function will be written as a delta function in the Bjorken variable x, which 
spikes at values of Xn corresponding to the excitation of the final state n. In the cm. system, 
this leads to the following exact transformation 

5{M^ - Po - go) = 5{x - ^n)-^ (4.24) 



with Xn the value of x at the bound state mass M„, given in Eq. ()4.26|1 below. 

In the DIS limit, the transverse currents dominate the scattering. The exact result for 
the transverse tensor, W^x-, that describes the scattering from the ground state with 7_ 
structure, to a single transverse final state n with 7-73; structure, is 

W..\=5{x-Xn)-^[j'^^-{Q\y) + Jl^-{Q\y)] , (4.25) 

where the reduced transverse current was defined in Eq. ()3.18|) . (Here the fi subscript is 
replaced by n, and denotes a transition between the ground and nth state.) Next we observe 
that the value of x at which the final state is excited is 

Q2 

= M2 - M2 + Q2 ' (4-26) 
and that the spacing between neighboring levels is 

A - - 

Q 

with = M'^j^^ — M^. Hence the exact result for the total inelastic tensor W^x is 



EH^-^,) A.,. (^) (^) {j-'i^.y)^J-'i^.y)}\ (4,28) 



\6Ml) \Xn+l 

where the sum is over all possible final states n. 
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2. Hadronic result in the DIS limit 



The DIS limit is defined by the requirement that Q"^ ^ oo with x held constant. In this 
limit both and become large, with 

Xn ^ M2%Q2 ' (as Q' and oo) . (4.29) 

Study of the DIS limit therefore requires estimating the transition form factors for large 
final state number n, but for all x (since Xn depends on Q^). In the DIS limit, y — > x [easily 
seen from the expansions ()E1|) ]. so y will be replaced by x in the following discussion. 

When n is large the "tails" of the wave functions are very small, and the exact transition 
current ()3.18|) can be approximated by neglecting the contribution from the region C, € [0, x] 
(where ^' < 0), giving 

l,x^^2 eimoiX /■! $„(e',P/)$(e,^) peff,„ ^2^ 



-■-»^i< ^"'^'1^riHi|'''' ^''-'^->' (4.30, 

where we continue to suppress the label n = 1 on the ground (initial) state wave function, 
and 

ri = -e (^) < (4.31) 

in the region of integration. When 77 < and Q"^ —>■ 00, the relation ()D10|) may be used to 
approximate the quark form factor by 

F^iv,Q')^ TwK] (4.32) 

TlQ^T] (1 - 1]) 

giving 

Ff{v,Q')^l. (4.33) 
Equation ()4.3()|) is now further reduced by expanding $(^, P)/^ around = x 

mn ^ Hx,p) ^ ^ _d_ / $(x,p) \ _^ 1 ^2 ( Hx,p) \ ^ _ _ _ ^434. 

^ X dx \ X J 2 dx^ \ X J 

where the higher order terms are proportional to powers of S 

6 = e'(l - x) . (4.35) 
This replacement gives the following series for J' 

$„(e',P/) 



Jn'^'iQ^x) = eimoi$(x,P) / dC 

Jo 



m=oo 



m=l 



^ e 

ml \dx I \ X 



+^.™o.EC^m (4.36) 
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where 

C= t di'r-'^n{i',Pf). (4.37) 

J 

Using the identity ()B10|) and the other estimates given in Appendix |Bl it follows that 

lim|:-.^-.0. (4.38) 

"^'^ Cn n 

Hence the large n approximation to ()4.36p is 

J■„l'-(Q^a;) ^ eimoi<l'(x,P) f 

Jo 4' 

= eiC„$(x,P). (4.39) 

The plus component of the transition current is smaller than the transverse current by a 
factor of rriQi/q^ [compare the two results in Eq. ()3.16|) ]. and hence is obtained immediately 

from the result ()4.39|) . Using Eq. ()E1|) for q_, and Eq. ()4.29|) to replace \J {1 — Xn)/xn, gives 

J:„^'+(Q^ x) ^ -V2e, Cn $(x, P) . (4.40) 

The minus component can also be obtained from Eq. ()4.39p . Reviewing the derivation, we 
see that the minus component, in the DIS limit, is obtained from the transverse component 
by replacing 



Hence, multiplying ()4.39|) by 

- v^e. C„ ( '"" + '-^J"'^' ) P) , (4-42) 

where the ^' integral was evaluated using identity ()B11|) . 

The transition currents can be shown to be gauge invariant by useing the Ward-Takahashi 
identity on the quark-photon vertex, and reducing the result using the Bethe-Salpeter equa- 
tion. 

These currents will now be used to determine the structure functions Wi and W2 in the 
hadron picture. 



3. The DIS cross section from the transition form factors 

The cross section in the DIS limit can now be obtained from the exact formula ()4.28p . 
As long as x is not too close to unity each term in the sum must correspond to some state 
with large n (the states with small n "pile up" near x = 1 where their total contribution 
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is very small), and the approximate result for the transverse current, ()4.39p . may be used. 
Furthermore, for large n, may be approximated by using ()B13|) . and Xn+i — Xn- Hence 
the final hadronic result [denoted / to distinguish it from the partonic result / of Eq. ()4.15p ] 
in the DIS limit is 

fix) = 2F,{x) = ^ = E '^(^ - ^n) |ei P) + (-l)"e2$(l -x,P)j . (4.43) 

Note that Fi{x) is a series of separate spikes, zero for all x except at particular values x„, 
where it is infinite. The partonic function, Fi{x), is a smooth function of x. To compare 
Fi{x) to Fi{x), we must average over x. This will lead to the concept of duality and will be 
discussed in the next subsection below. 

Now look at the component of the tensor. Replacing the transverse components of 

the current in Eq. ()4.28|) with plus and minus components, and taking the DIS limit, gives 

^oo n 

X {j^''{Q\x) + J^'-{Q\x)} {J-„^'+(Q^x) + J-„^'+(Q^x)}* . (4.44) 

Substituting Eqs. ()4.4U|) and ()4.42p . and droping the term proportional to (—1)"^^ (which 
averages to zero) gives 

W^-+le,te.n. = E " ^n) 6? ^^(x, P) . (4.45) 

V n 

This is a nonleading term and can be neglected in the DIS limit. Hence the hadronic picture 
also gives the Callan-Gross relation, and it is sufficient to compare the functions Fi and F\ 
only. 

If Wi were zero, Eq. ()4.45|) would give the following scaling relation for v'^W2 



n2 ™2 

^'^2L,,erm - ^^-+le.term = 2^ /(^) le, term • (4-46) 



This is to be compared with ()4.22|) . Hence, the duality of Wi and the duality of the non- 
leading terms in v'^W2 depends on comparison of the same functions, f{x) and f{x). 
We discuss this comparison now. 



D. Duality and its implications 

In any physical measurement, the detectors which define the final state will accept a range 
of final electron energies 5E' . For fixed Q^, this can be converted into the acceptance of a 
range of values of x centered at Xi with width 5x so that x varies over the interval bounded 
by Xi± = Xi 5x/2. Then the experimentally measured hadronic structure function in the 
DIS limit can be computed from ()4.43|) 



^' 5x 



dxf{x) 



— y Ax 
Sx ^ 



ei$(x„, P) + {-ir-'e2 ^l-Xn,P) 



(4.47) 
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M 

average 

FIG. 12: Quark- hadron duality 



Note that this is the sum of smooth terms proportional to ef and e|, and a rapidly oscillating 
interference term proportional to 6162- In the DIS limit, the separation A^. between states 
approaches zero, and therefore a large number of states are necessarily included in any 
interval 6x. The interference term, which changes sign as (—1)""^, therefore averages to 
zero. This cancellation of the interference term was pointed out originally by Einhorn 
and empahsized recently by Close and Isgur in the context of the nonrelativistic quark model 
0. 

The smooth terms can be approximated by their value at the center of the interval. Using 
the fact that the number of states in the interval times the x spacing between them must 
necessarily equal the width 6x, so that 

E = ■ (4-48) 

we can reduce ()4.47|) to 

(4^ ^ {e?$2(x„P)+e^$2(l-a:„P)}l ^ A.„ 

+2eie2 ^ E {-ir~'KM^n, P)Hl - Xn, P) 

^X n€Xi 

^ e?$2(x„ P) + el<!>\l - x„ P) = f{x,) , (4.49) 

where f{x) was defined in Eq. (j4.15|) . 

This result should be compared with the partonic result, averaged in the same way. Since 
the partonic result is smooth (and the interference term already neglected) the partonic 
average is trival 

(/) = ^ n dxfix) ~ f{xi) ^ n dx = f{xi) . (4.50) 

^ ' dx Jxi- OX Jxi- 

Hence we see that 
proving duality. 
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Note that a similar identity holds for the nth moments of /. In the high limit when 
the final states are numerous and close enough together so that A^.^ (ix — > 0, the mth 
moment of the hadronic structure function is 



Jo 



el <l>2(x„,P) + e^$2(l-x„,P) 
el $2(x,P) + $2(1 -x,P) 



dx x'^fix) 



(4.52) 



Quark-hadron duality means that properly averaged hadronic observables in a certain 
kinematic regime (high Q"^) can be described by perturbative QCD, as schematically repre- 
sented in Fig. ^1 and demonstrated in Eq. ()4.51|) . Recent data from Jefferson Laboratory 
exhibit duality, and this has focused attention on this old subject. Recent studies show that 



luDject. m 

the weak decays within the 't Hooft model exhibit duality 0| ^ and in the framework 
of DIS, Close and Isgur jl^ studied its emergence in the context of realistic nonrelativis- 
tic quark models. A recent paper by Isgur, Jeschonnek, Melnitchouk and Van Orden 
studied duality using a relativistic Klein-Gordon equation with a confining interaction. This 
work is complementary to these other studies. 

While we have proved that duality emerges in the DIS limit, and have corrected the 
pathologies of earlier work, we have not yet studied how duality emerges as a function of 
Q^. A future paper will present numerical studies of the onset of duality for rQCD ,15i] . 



V. CONCLUSIONS 



We have shown that a very satisfactory description of deep inelastic scattering can be 
obtained from large Nc QCD in 1+1 dimensions, provided the quark spin degrees of freedom 
are extended to the full 1+3 dimensions (along with the leptons and the electromagnetic 
field). We refer to this new model as reduced QCD, and suggest that rQCD may be used 
to provide insight into the behavior of many other physical processes. In the applications 
to DIS discussed in this paper, we are able to confirm the Callan-Gross relation, and that 
duality emerges in the deep inelastic limit. 

Our theoretical discussion teaches several lessons. Electromagnetic gauge invariance and 
duality emerge only if the quark current is fully and consistently dressed. The way in which a 
quark form factor ernerges from this dynamical model mimics the phenomenlogical approach 
of Gross and Riska [iSl, showing that their prescription has a dynamical justification. Our 
treatment of confinement removes all singularities, giving a finite mass for the dressed quark. 
This mass is gauge dependent, showing that it is not a physical quantity, which can be taken 
as a demonstration of confinement. But the fact that we have finte dressed quark masses 
in the presence of confinement provides justification for the constituent quark model. The 
simultaneous presence of confinement and finite quark mass is possible only if the "on-shell" 
quark scattering amplitude is identicallyt zero, and we have shown that this is indeed true, 
as modeled previously in Ref. [Tt^ . 

The treatment of reduced QCD presented here lays a foundation for much further study. 
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APPENDIX A: REDUCTION OF THE TWO QUARK INTERACTION 

Integrals over the variable /c+ like the one encountered in Eq. ()2.35|) appear several places 
in this paper. Here we evaluate it in full generality. 
Consider the integral 



Substituting for V gives 



d^k V{p,k)k-(k- -r_)F(k_,r) 
(27r)2 [mf — k"^ — it] [m^ — {k — r)^ — ie] 



(Al) 



dPk 



F{k_,r) 



F{p_,r) 



TT^ J {p_ - k_y [d2{k-,k+,r) d2{p-,k+,r) 
where the denominator d2 is 



(A2) 



d2{k_,k+,r) 



m 



i -2^4 



— z 



k. 



[k — r)- 



m 



yr- 



^-2k, 



— I 



yr- 



(1 -y) r_ 



2{k-r)+-i 



2{k - r) 



[k — r)- 



[l-y) r_ 



(A3) 



and we have introduced the momentum fractions ()2.3(i|l . This shows that the denominator 
has only two poles in and that they will both be in the same half of the complex plane 
(giving zero for the integral) unless 



< 1/ < 1 . 



(A4) 



Doing the integral over and expressing the answer in terms of the momentum fractions 
y and z gives 



l2{zr_,r) l2{z,r) 



9^ /■! dy j F{yr_,r) 
n Jo {z- yf \ A(?/,r2) 

-e{z)e{i-z/^^''-''^ 



:z)e{i-z) ^^''-''^ 



TT A(z, r^l 



roo 

+ / 

-oo Jl 



A(z,r2) 
dy 



(z-yy 



(A5) 



where the new denominator A{y, r^) is 



A{y,r^ 



+ 



rrio 



y (l-y) 



(A6) 



31 



[Note that this result holds for both signs of r_. This is because the sign of the countor 
integral, which depends on the sign of r_, is cancelled by the sign of the integral over 
/c_ = ?/r_, which must be changed from J^* — when r„ < 0.] It is important to 
note that the integral over y is not restricted in the term in ()A5|) that is proportional to 
F{zr_,r) because for that term the restriction ()A4|) applies to z and not y. This term 
not only regulates the singularity in the first term (which arises only if < z < 1), but 
also contibutes an additional contribution. Furthermore, if z lies outside of the interval 
< 2; < 1, the first term is non singular (except at the end points of the interval) and is 
nonzero! Hence, in the region < 2; < 1 the integral ()A5|) becomes 



dy i F(yr^,r) F(zr^,r)\ g'^ Fizr^ 

h[z,r) = — ' 



TT Jo {z — yY \ A(?/,r2) A(2;, r^) j -n z{l — z) ^{z^r"^) 

\T( \F{yr^,r) g^jix F{zr_,r) ^ m m 

A[y,r'^) z[l — z) A(z, r^j 

and outside of this region it is 

T( \ 9^ r F{yr-,r) F{yr^,r^) ^ rn m f \q\ 



where we have defined 



vr 



^2 



(A9) 



y'f 

The multiplication by in ()A9|) is needed to scale the momenta from k_ y, etc. 

We emphasize that the identities ()A7|) and ()A8|) hold for both r_ > and r_ < 0. The 
first case is needed for the bound state equation and the second for electron scattering. 



APPENDIX B: IDENTITIES INVOLVING BOUND STATE WAVE FUNCTIONS 

In this appendix we derive a number of identities that hold when $^(0, r) = $„(!, r) = 
(i.e. in all cases but the chiral limit). Equation (j2.42|) leads immediately to the following 
identity 



^ m, ^ j ^^^^^ ^ ^ _ 

\ z I — z 



Integrating by parts gives 

/o ^ Jo il-z 



mil ^^^"^ = ^1 d^TTZ^- (B2) 



An additional identity is derived by multiplying ()2.42|) by 2 (1 — z) ^n{z, r) and integrating. 
First observe that 



2 ['dz{l-z)^'^{z,r)-f'dy- 
Jo Jo [z — y) 

= /' /' dz dy $;(^, r) r) ^ = . (B3) 

JO JO z — V 
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Hence, integrating by parts and using the normalization condition ()2.40|) for the wave func- 



tion, we obtain 



Similarly we can prove that 



[l-z) 



(B4) 



(B5) 



Additional identities are needed for the study of duality. First, integrate Eq. ()2.37|) over 
z to obtain 



m, 



01 



+ 



Ml I dz^n{z,r) = [ dz 
Jo Jo 

Next, the commutator of the operator 

/C$„(x,r)=:f dy 

Jo 



02 



1 - z 



(B6) 



y-x 



(B7) 



with the Hamiltonian defined in ()2.37p is 



[H,]C] Mx,r) 



dy 



m, 



01 



m. 



02 



xy {l-x){l-y) 



(B8) 



as discussed in Refs. j2[ 0|. Multiplying this by r) and integtating over x gives zero 

on the l.h.s., and hence the relation 



01 



m. 



02 



1 - z 



(B9) 



Taking the square root of both sides of this equation gives 



moi / dz = ±mo2 / dz— = 6„ , 

Jo z Jo L — z 



(BIO) 



where the sign is positive if the phase of $„(0, r) is the same as $n(l, r) and negative if it is 
opposite. Since the states are non-degenerate, the phase of the nth eigenstate is (— l)*^""^). 
Combining this with Eq. ()B6|) gives 



Ml 



dz ^n{z,r) 



moi 



"^02 



Cn 



(Bll) 



where the constant C„ is yet to be determined. In this paper we have chosen the sign of 
$^(2;, r) so that it is always positive as z — > 0, so that C„ > for all n. 
Our demonstration of duality depends upon the relations 



hm Cn — > C, 

n—>co 



00 ^ gv^ 

\\m (mI^, - Ml) -.Cl = T^g" 



(B12) 
(B13) 
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which can be derived approximately using arguments given by 't Hooft and Callan, 
Coote, and Gross [2], given here for completeness. First, for large n, the boundary conditions 
require that the normalized wave function go like 

sin{mTz) . (B14) 

One test of this approximation is shown in Fig. E| If n is very large, so that M„ 3> rriQi, the 
behavior of the wave function at the endpoints can be ignored, and near z = 1/2 Eq. (|2.42p 
can be approximated 



M^^r.{z,r)--^4- (B15) 
TT j-00 y — z 



Substituting the ansatz ()B14|) into this equation gives 

, ^9 . , , ■) cosfnvr-u) 9 /■°° cosinny) 

Ml sMmiz) = -g^nf ^ — — = -g^nRe / ^ — — 

J -00 y — z J -00 y — z — le 



= ng n sin(7r?7,z) . (B16) 

Hence — >■ irg'^n as n — >• cxo and identity ()B13|) is proved. To prove ()B12|) we consider the 
matrix element of (IB8I) 



Jo 

= (Ml, - M^) f dxdy^^^±^^^^2^^^^ . (B17) 

^ ^ JO y — X 

where we used the fact that the phases of and $n+i are opposite. The integral is evaluated 
in the large n approximation 

/ d^dy^n^^M^n{y.r) ^ 2 r ^^sin([n + l]vrx)sin(n.y) 

JO y — X Jo J -00 y — X — ie 

= 2 c/x sin([n + llvrx) cos(n7rx) 
Jo 



Hence, 



^ ng^ ~ Cl (B19) 

and ()B12|) is proved. 

A numerical demonstration of the identities ()B10|) . ()B12|) . and ()B13|) is given in Figs. IT^ 
The identities are satisfied to better than a few percent for state numbers between 40 and 
80, provided the number of splines is at least 160 (twice the number of the maximum state 
of interest). 

Two additional relations that follow from the completeness relation ()2.5H) and the bound 
state equation are needed. U z E [0, 1], then: 

J2 ^niz, r„) /' dz' $„(^', r„) = 1 (B20) 
Jo 

E f ^ + 1^ - M^] '^niz, r„) dz' r„) = . (B21) 

n \ Z l-Z ) Jo 
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FIG. 13: Numerical test of the relations IjBlOp . (|B12|) . and IIB13|l for moi = 1.5 and mo2 — 3 and different 
numbers of splines, Us- Solid circles are values of the l.h.s. of (|b10)| divided by Ctx^; open circles are the 
absolute value of the r.h.s. divided by Coo- They should be equal for each state, and approach unity as 
n ^ oo. The tiny solid squares are the values of {fJ-n+i — Mn)/""^ which, according (|B13|) . approach unity as 
71 — > oo. The open triangles are the square of the normalization constant divided by 2 (for wave functions 
initially normalized with their maximum value equal to unity) , and approaches unity if the approximation 
(IB 141) is valid. 



The proof of the second relation ()B21|) follows from application of the bound state equation 



E 



I — z J Jo 



EL dy 

n 

E 



+ 1 +1 My 



z — e J z-\-e 



{y - zf 

YX\^ dyKiz,rn) I dz'^n{z',rn) = 

„ Z Jz-e Jo 

:5]$::(z,r„) ['dz'^^{z',rr,)^0, 
„ Jo 



dz'^n{z\rn) 



(B22) 
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where the identity ()B20|) is used to ehminate the nonsingular integrals over the intervals 
[0,z — e] and [z + e, 1], and the singular part of the integral is proportional to e which 
vanishes as e ^ 0. 

If Eq. (|2.4ip is used, modifications of the completeness relation (j2.5H) and the identity 
()B20|) . which hold in the region z ^ [0, 1], can be derived. The results are 

E f- + - M^] ^niz, r„) /' dz' r„) = • (B24) 

n \ z 1 — z J Jo z{l — z) 

Note that the completeness relation and all of these identities hold for any two quark 
Hamiltonian. In particular, the Hamiltonian may describe bound states with flavor (where, 
in general, mi 7^ as well as flavorless bound states where mi = 1712. In some applications 
we will study the coupling of a photon to a qq pair. Such a coupling necessarily involves 
flavorless states only, and we will use the identities for the flavorless states. 



APPENDIX C: SUMMATION OF THE QUARK SCATTERING SERIES 

Here we show how to sum this series ()2.58j) and obtain a useful form for M. 

Since V is independent of the plus components of momentum, the integrations over the 
plus components of momentum can be carried out. Using the identity ()A7|) [with g_ < in 
this case], the first term on the rhs of ()2.58p becomes 



{MGO),ip';q) = I 



(fp V{p',p) 0{p^, q) 
7r2 d2{p-,p+,q) 

0{p.,q) 0{p'_,q) 



L ^^^^"^ ' A(2, g2) vr z'{l - z') A(^', q^) 
Jo \ ^ ^ Tizil - z) I A(z,q^) ^ ^ ^ ^ 



or 



(MGO).(p';,) ^ . / fp yiP'fO(P-f , /',,Vo(/..)g|g::l|l if ^ |0. 1] . (C2) 

J TT^ d2[p-,p+,q) Jo A{z,q^) 

Here the replacement ()A9|) has been used to convert V ^ V (or Vq). Similarly, the second 
term becomes 

d^k V{p',k) /■! . gH{y- z)] 0{p^,q) 



tt'^ d2{k^,k^,q) Jo [ ' n z{l — z) j A{z,q 



■2\ 



7ri/(l -y) j A{y,q^) 

, dz I V y, z) \ — ^ liz e 0, 1 , C3 

I -K z[l- z) \ A[z,q^) 
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or 



{MGO)^{p'-q) = l'dyVo{z',y) 



r^4v(.,.)-^^^|^ If .'^[0,1], (C4) 
Jo [ n z{l — z) ) A[z, q^) 



and the third term follows the same pattern. Replacing {q^)'^M{p' ,p; q)/^ M{z' , z; q"^) 
and 0(2;, g^), the series fj2.58|) . in the interval z' G [0, 1], is clearly summed by 

the following equation 

^ A/TM 2^0(^^ , , g'6{z'-z)]0{z,q') 

dzM[z,z;q)— —= / dz<\{z,z) ' 



^ ' ' A{z,q^) Jo I ' ' ' nz{l-z) jA{z,q^) 

Miz',y;q') j g'Sjy-z) ] 0{z,q') 

Jo A{y,q^) Jo [ n z[l — z) J A{z,q'') 

Introducing the shorthand notation 

-(^'•-^)-^. («) 

Eq. ()C5|) can be rearranged into the following form 



f\zF{z',z;q^)\lA{z,q'')+ 

Jo [\ z{l-z)J A{z,q^) Jo A{y,q 

= f dzF{z',z-q') (H{z)-q') ^ 

Jo ^ ^ A[z,q 



I -Kzil-z) ] A{z,q^) 



which requires that F satisfy the equation 

F{z', z, q') {h{z) - q') = V{z', z) - ll^l^ . (C8) 

^ ' 71 Z[i — Zj 

This equation is easily solved with the use of the Greens function: 

F{z\z-q') = -[dy ^) " T^|r^} ^'(y, ^, g^) 

\y My,r„)$„(.,r„) 



Putting it all together gives Eq. ()2.59|1 . 



i ' 7r?/(l-?/) J ^ g^-M2 



g2 - M2 



(C9) 



37 



If -z' ^ [0, 1], the series reduces to the same formal result 



dzM{z',z;q^ 

Jo 



2^ 



A{z,q'^) Jo ' A{z,q 

+ ['dyVoiz',y)—^ f'dzMiy,z;q'^ 



A{y,q^)Jo ' A{z,q') 



f dzVo{z',z) ^ f dyVo{z',y) 

Jo A[z,q'^) „ Jo 

dz^n{z,rn) 0{z,q 



A{z,q^ 



= E A{z\ Ml) /' dz $„(^, r„) 0(^, g^) , (CIO) 

where the completeness relation ()2.5H) and Eq. ()2.4H1 were used in the last step. This shows 
that (j2.59p can be applied for any z' G [—00,00]. 



APPENDIX D: SUMMATION OF THE QUARK CURRENT SERIES 



The series flH.7|l for the dressed quark current is evaluated in this Appendix. The "first" 
integral over /c+ in each term in the series is evaluated as in Eq. ()C1|) or ()C2|) 



d'^k 
(2^ 



V{p\k)-i^S,{k)-i^S.,{k-q)-i_=i 



d^k V{p',k)Nt{k-,q-) 
7r2 dl{k_,k+,q) 



dy ^Viz',y) 

dy yoiz',y) 



g^6{z'-y) \ Nt{k^,q^ 



?/(!-?/) j A,(?/,g2) 
Nt{k.,q.) 



where 



iVf(fc_,g_) 



7_[mo^ + fc_7+] Y [mpi + (fc_ - q-h+] 7 
4k_{k_ - qJ) 



if z' E [0, 1] 
if z' i [0, 1] 



(Dl) 



:d2) 



and d\ and Aj are ()A3|1 and ()A6j) with equal masses m\ = m2 = rrii. Note that the only 
dependence of the integrand on A;+ was in the denominator d\, so the /c+ integral could be 
evaluated using the methods of Appendix O The second integral is similarily evaluated. If 
z' G [0, 1], the following result is obtained for the series Eq. ()3.7|) 



Jt{p\p' -q)=jt{z\q')=r+ dy V(/,y) 



+ / dy'l V(z', y' 



gH{z'-y)\ Nf(y,g) 



7rz/(l - y) J Ai{y\q'^) 



9'6{z'-y') \ 1 
ny'il-y') J A,(|/',g2) 7o 



dy \v{y',y) 



+ 



7^^+ / dzM{z',y;q'' 



Kiy,q) 
My,Q') 



gH{y'-y) \ Nr(y,g) 
(D3) 
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where Nf(|/,g) = N^ik^jqJ), and the series is summed using Eq. ()C5|) for the scattering 
matrix M. Using the solution ()2.59|) for the scattering matrix the result, for all values of 
z' = p'_lq_^ is 

jr(^', q') = 7^ + E Ml) %^ /' dz r„) Nf (z, q) , (D4) 

where ji{p', p'—q) = jti^', q^)-, and Aj is the A of Eq. ()A6|) with equal masses nii = 1112 = rrii. 

The same expression also holds for z' ^ [0, 1], but the derivation differs. As before, the 
full result can be deduced from the form of the first two terms of the series (|3.7|) . Guided 
by ()D1|) . and using the steps already displayed in Eq. ()C10|) . we obtain 



A,(|/',g2 

+ / dyYQ[z,y)-—^— / dy \Y[y ,y) ^\^, ^ 

Jo I\Av\a^) Jo TTvil-v) \ AAv.a^] 



/^i{y',q^) Jo \ ' 7r?/(l - ?/) J Ai(?/,g2) 

Nr(i/,g) 



+ ■■■ =Y+ I dzM{z',y;q 
Jo 

= 7^ + EA,(;.',M^)M^ /^^<|.„(^,r„)Nf(^,g). (D5) 
n M^-Q Jo 

We see, as we did with Eq. flTHHll . that (IDill holds for all z' . 

Equation ()D4|) expresses the dressed current as a sum over contributions from the bound 
states that couple to the photon, showing that vector dominance is a rigorous consequence 
of this model. 

The structure of the numerator Nf depends on the bare current: 



7- : 


Ni (y, g) = 7- 


7+ : 


N+(i/,g)= -7-- 


Ix ■■ 


Nf(2/,g) = -7x7- 



2g2^(l-y) 
" 2q_y{\-y) ' 



(D6) 



Note the interesting structure of the correction term to the bare current 7^;, which will be 
of central importance in our discussion below. Substituting these results into ()D4|1 gives the 
following 



{ n q JO ) 

We organize these expressions by introducing the quark form factor, defined in Eq. ()3.8p . 
This immediately gives Eq. ()3.9p for the j_ current. The other currents in ()D7|) can also 
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be written in terms of Fj. To this end note that the bound state equation for equal mass 
quarks imphes that 



Ml / rfy r„) = m2. / dy 



/o y{l-y) ' 

Hence the quark mass can be removed from the second term in j+, giving 



(D8) 



9i{z',q) 



m, 



2qi 



dy- 



Ml -q^ Jo y{l - y) J 



q- q^ 



5:A,(/,M2)$„(/,r„) f^dy^n{y,rA , (D9) 



where gi{z',q) defined by this equation should not be confused with the two body Greens 
function or the bound state vertex function. If z' G [0, 1], the second term can be reduced 
using the identities (P2n|l and (P2T|l 



E 



m. 



Oi 



z'(l - z' 



z'(l - z' 



^n{z',rn) / dy^n{y,rr. 



z'{l - z') ■ 



(DIO) 



The same result is obtained directly from identity ()B24j) if z' ^ [0, 1]. Hence, for all z' 

-2 1 „ , , g' 1 



g,{z',q) = -^F,(z',Q^) + 



' ^ i. 

27rg2 z'{l- z') 



~^-±F,iz',Q')-^^. (Dll) 



Finally, the component of the current can be similarly reduced. The second term in 
jx is proportional to 



Hiiz',q) 



M2-g2 J, yy(l-y) 



2g_ nioi 



J2MlMz',Ml 



dy^niy,rr, 



2q_ mo- 



Ml + Q^ Jo 

-|^A,(;2',M„2)$„(;2',r„) j\y^^{y,rA. (D12) 



The part of the second term in {} is the same constant derived above, and combining the 
results gives Eq. 1)3.111) . 



APPENDIX E: KINEMATICS FOR DIS 



Kinematics for DIS in the cm. frame are given in this Appendix. When is very large, 
the components of P and q in this frame are (see Ref. [l^ ) 



X 



2x 



\ Mlx{l-2x) 
2Q\l-x) 



40 



q- 



-Q^ 



X 



2(1 -x) 



x] 



Q V 2 
Q 

2x(l -x) 



M^x(l - 2x) 



2g2(l_a;) 

M2a;(l - 2x2) 
2Q^{l-x) 



+ ... 



1 + 



M2x(1 - 2x2 

2g2(i-x) 



+ 



(El) 



The energy and momentum of particles 1 and 2 in the final state are 



Pz 



Pw 



P20 



± 



Q 1 



X 



^ ^ x(M2 - 2ml - 2mi) ^ 



mf + p2 = — 



m2 + = — - 



2Q2(l-a;) 




x(M2 + 2ml - 2^ 



2Q2(i-x) 



+ 



x(M2 - 2m2 + 2m2) 
^+ 2Q2(i_a;) + 



(E2) 



There are two possibilities corresponding to the two terms in the current (j4.1(J|) . The first 
term in the current will be large only if the momentum of particle 1 is in the direction of 
(i.e. piz = \pz\ > 0). [Choosing it in the opposite direction requires a very large momentum 
fiow through the wave function, and is suppressed.] In this case 



Pi+ 



p'^1 



P2+ 



P^l 




x(M2-2mi) 
2Q^{l-x) 

x{M^ - 2m. 



2g2(i 



X 



x(M2 



2g2(i 

2mj] 



2Q2(i-x) 



+ ■■ 



+ 



x{M^-2ml) ^ 



(E3) 



where the superscript (1) is a reminder that these relations hold only for the first term 
in the current. The momentum relations for the second term in the current ()4.10|) follow 
by interchanging 1 and 2 on both sides of these equations. We see that the first term in 
the current gives pi+ Q and p2+ — > 1/Q, while the second term gives pi+ — > 1/Q and 
P2+ — >■ Q- The two terms describe kinematically distinct regions of phase space, and their 
interference can safely be neglected when the total cross section is computed. Using the 
definitions (j4.7p and (|4.9p . z and z' can be easily related to the Bjorken variable x. We have 



z = 1 



z' = 1- 



pjl 

p_ 

(2) 

Pl- 



1 — (1 — x) = X 



1 — (1 — x) = X - 



(E4) 
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APPENDIX F: THE BOUND STATE TRANSITION CURRENTS 



This Appendix includes details of the evaluation of the transition currents shown in Fig. 
[TUl The Ci term will be evaluated first. The 62 term can then be obtained by a simple 
substitution, as discussed at the end of this Appendix. 

The Ci term, given in Eq. (j3.13|) . is reduced by introducing the momentum fractions 



(P + A;)_ 



V 



so that 



f- 



P- q- ' Pf- i-y 

In this notation the x-type transverse quark current, evaluated in Eq. (jH.llll . is 

,2 



(Fl) 



(F2) 



Ix - Ix 1- 



9 



q- 



27rmoi (P/_ + fc_)(P_ + A;_) 
-fx + 7x7-Hi{ri,q) , 



q+ 

moi 



(F3) 



First evaluate the numerators, iV^, of the traces in ()3.13|) . The numerator is 



Tr 



l-lx{moi + {Pf + A;)-7+} ilx + Ix 7- ?)) 

X {moi + (P + A;)_7+}7_ (/i;-7+) 
16 A;_ [moi g_ + 2(Pj + k)^{P + k)_Hi{r], q)] 



16 k-Q- 



q+ 



32{Pf + k)4P + k)^^Fi{7], Q 



moi moi 



moi 



(F4) 



where the trace has been evaluated in a Dirac space of four dimensions. The numerator of 
the minus component of the current has only one term 



Tr 



l-{Pf + ^)"7+7-(P + ^)-7+7- (^-7+) 



32 k.{Pf + k).{P + A;)„ [1 + Pi(r7, 



(F5) 



The trace for the plus component is 



Tr 



7_{moi + {Pf + A;) -7+} (7+ + 7_ Gi{ri, q)) 
X {moi + (P + A;)_7+}7_ (A;_7+) 
16A;_m^ - 32A;_(P/ + A;)_(P + A;)_ ^ Fi{r], Q^) 



IQk^ [mi - 7]{1 - T]) Q'Fi{r],Q' 



(F6) 
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Inserting the numerator ()F5|) into the general result ()3.13p gives the following result for 
the minus component of the current 



{f.\J4Pf,P)\t.) 



ei term 



32iei 
32zei 



d^k {Pf + A;)_(P + k).k_ Gf {-Pj - k, -k) d (P + k, k) 
(27r)2 [ml - {Pf + A;)2 - ie] \m\ - {P + kf - le] [ml ~ k^ - le] 

"^'^ ^ Gf {-Pf - k, -k) Gi {P + k, k) \l + Fi(77, 



(271)2 ^3 



(F7) 



where the subscript - on /_ reminds us that the final state must have a pure 7_ structure (if 
the ground state does), is the product of the denominators of the three quark propagators, 
with three poles in fc+ 



d-i 



f- 



/- 



Mf 

-j^-2k+- let 



(C-1)^- 



— 2fc I — ier 



(F8) 



and the e's change sign according to 



e. = e/{Pf + k)^ = e/{ePf- 
e, = e/(P + fc)_ = e/(eP-) 
ec = 



^/[i^-y)P- 



eA„ = e/[(e-l)P-]. 



(F9) 



Since the vertex functions G do not depend on k+, we can evaluate the k+ integral. It will 
be nonzero only when the three poles of ()F8|1 do not all lie in the same half plane. Since 
g_ < for electron scattering, P > Pf_ > 0, which implies that 



Hence 



<y <1 



e-i<e-y<e: 



(FIO) 



(Fll) 



and all three poles will be in the same half of the complex plane (giving zero for the integral) 
unless 



< ^ < 1. 



(F12) 



There are two terms, depending on the sign of C, — y- Closing the k^ contour in the upper 
half plane, and using dk_ = P_ dC, and {C, — 1)P- = — 1)-?/- gives 



{f_\J4Pf,P)\^_) 



ei 

term 



16ei P^Pf 



1 d^ Gf{C,Pf)Gd^,P) 
(27r) A(e,p2)A(e,P2) 



Wo (27r 



d^ Gf{C,Pf) G.(^,P)[l + Pi(r^,g2)] 



A(e,P|) (A(e,p2) 



A(e,p|: 



, (F13) 



where the first term is the contribution from the pole due to the zero in the third term 
in ()F8|1 above, and the second from the pole due to the zero in the first term. Using the 
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definition ()2.39p of the wave function and combining the two contributions in the region 
[0, y] we get Eq. TTT^ with 



7^ 



A(e,p?)-(i-y)A(^,p2) 





e'(i-m|-(i-eK 






(1-e) {M]-{l-y)Mf) 





if M} 



oo . 



(F14) 



This factor is needed for positive values of in the interval [0,?/] and negative values of ^' 
in the interval [—y/{l — ?/),0]. In this region the denominator has a zero only if m\ < 0. 
In order to avoid the discussion of such cases we limit numerical applications to cases with 



m 



2>0. 



Results for the other components of the transition current, Eq. are obtained by a 

similar arguement using ()F4|) and ()F6|) in place of ()F5|) . 

Now consider the modifications required in order to evaluate the 62 term. Using the 
momenta defined in Fig.^[ and the same definitions of momentum fractions ()F1|) and ()F2|) . 
the 62 term is obtained from the ei term simply by substituting mi ^ m2 and ei 62- 
However, the momentum fractions in the wave functions are, by convention, the fraction 
of the momentum carried by the quark mi, and the momentum fraction in the quark form 
factor is that of the the outgoing quark, and hence these functions must be written in terms 
of 



6 

V2 



-k. 



p. 



f- 



-jP + k). 

q- 



(F15) 



These observations lead immediately to the final result ()3.17|) . 



APPENDIX G: THE MODIFIED CUBIC SPLINES 

The two body equations were solved using a modification of the standard cubic splines 
employed previously in many problems. The standard splines are defined on 4 segments of 
length h, bounded by the 5 points a = (n — 2)h, b = a + h,c = b + h,d = c + h,e = d + h: 



1 

— < 

4 



(x—a)^ 
1+3 

1 + 3 

(e— a::)'^ 
/i3 



(x-b) I (x-b)^ 
h ~^ h-2 

jd-x) I {d~xf 
h ~^ h-2 



{x-bf 



{d-xf 



if a < X < b 
if b < X < c 
if c < X < d 
if d < X < e 



(Gl) 



If the interval [0,1] is spanned by — 2 standard splines, there must be 
length h = l/{ns + 1). The first spline is numbered n = 2 beginning at x 
is numbered — 1 ending at x = 1. 



1 segments, of 
and the last 
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2.0 



2.0 



i'Wi.O - 




Y{x) 1.0 - 




0.5 



0.0 



FIG. 14: The modified cubic splines for the minimal case of = 3, with four segments. In the left figure 
the left- and right-most splines at the boundary approach as a:°'^™^ and (1 — a;)*^-^^"^ appropriate to quark 
masses moi = mo2 = 1-5. In the right figure the bare quark masses are zero and the splines approach a 
constant at the boundaries. Note that the sum of the splines (the dotted line) is a constant in the central 
region, and throughout the whole region in the chiral case. 



The modified splines used in this paper consist of the standard splines plus one additional, 
non-standard spline inserted at the beginning and end of the interval [0,1]. These non- 
standard splines are defined on only three segments, and will be numbered 1 and Ug- They 
are 



Si{x) 



{2h~x) 
h 



1 + 3 



+ 



3/? 

{2h-xf 

h2 



{2h-xf 
h3 



(■ih-xf 
4h3 

(x-l+'ihf 



di 



4^3 




^ 4 


'(l~h~x) 


h 


{1-x 
L 1,0. 





{1-h-x)^ 
. / 



{l-h-xf 



/l2 



h3 



if < x < /i 
li h < X <2h 
if2h<x<3h 

iil-3h<x <l-2h 
ifl-2h<x<l-h 
iil-h<x<l . 



(G2) 



where the (3i are the fractional exponents determined by the boundary conditions [given in 
Eq. ()2.49|) ]. and the coefficients di depend on the exponents (3 



di 



(3-/3)(2-/?)' 



d. 



3(1-/3) 
(2-/3) 



(2/? - 3) 
(3-/3) 



(G3) 



where (3 = Pi for or /?2 for Sn^. These coefficients were fixed by the requirement that 
the spline and its first two derivatives be continuous. Figure IT^ shows two examples of the 
modified splines for the minimal number = 3. In this case the central spline spans the 
full interval [0,1] and has 4 segments. 

Using these splines as a basis, the equation is reduced to a matrix equation that is solved 
with the standard eigenvalue subroutine packages. 
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